Noiseless Collective Motion out of Noisy Chaos 



Tatsuo Shibata 1 , Tsuyoshi Chawanya 2 , and Kunihiko Kaneko 1 
1 Department of Pure and Applied Sciences, University of Tokyo, Komaba, Meguro-ku, Tokyo 

153-8902, Japan 

2 Department of Physics, Graduate School of Science, Osaka University, Toyonaka 560-0043, 

Japan 

(November 25, 1998; submitted to Physical Review Letters) 



Abstract 



We consider the effect of microscopic external noise on the collective motion 
of a globally coupled map in fully desynchronized states. Without the exter- 
nal noise a macroscopic variable shows high-dimensional chaos distinguishable 
from random motions. With the increase of external noise intensity, the col- 
lective motion is successively simplified. The number of effective degrees of 
freedom in the collective motion is found to decrease as — log a 2 with the ex- 
ternal noise variance a 2 . It is shown how the microscopic noise can suppress 
the number of degrees of freedom at a macroscopic level. 
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Chaotic motions have been observed experimentally in physical, chemical and biologi- 
cal systems. Since the evolution of these systems is subjected to external fluctuations, the 
observability of deterministic chaos depends on how the external fluctuations influence on 
it 0,0- Motivated by this point, extensive studies have been carried out about the enhance- 
ment of predictability and unpredictability of chaotic motions |fJQ . 

So far such studies are restricted to low-dimensional dynamical systems. Low- 
dimensional chaotic motion often arises as a macroscopic motion out of microscopic chaos 
with many degrees of freedom. Let us consider external fluctuations imposed on the mi- 
croscopic level rather than the macroscopic level, as is probable in natural systems, such as 
fluid turbulence, or neural systems with a large number of neurons. Since chaos can am- 
plify a small-scale error, it would be natural to ask a question how such a low-dimensional 
macroscopic chaos is possible out of high-dimensional chaotic systems subjected by external 
fluctuations. 

To address the question, we note that in certain coupled dynamical systems a macroscopic 
variable shows seemingly low-dimensional motions, while microscopic variables keep high 
dimensional chaos. Such a phenomenon has been extensively studied as a collective motion 
in coupled map lattice ||, globally coupled oscillators ||, and globally coupled map |6|-|14|. 
In the present Letter, we focus on the effect of noise on the collective motion of globally 
coupled map (GCM). 

The present GCM consists of N elements iterated by a local dynamics f(x) with a global 
coupling among elements and an external noise. The dynamics is given by 

e N 

x n+1 (i) = (1 - e)f(x n (i)) + fMj)) + £«(*), (1) 

iV 3=1 

for the z'th element at time step n. Here, we adopt the logistic map f(x) = l — ax 2 for the 
local dynamics, and Gaussian random process for £ n (i), with = and (£ n (i)£m(J)) = 

a 2 S nm Sij. The variance of Gaussian distribution is denoted by cr 2 [[TJJ. 

In the noiseless case (a = 0), if the coupling strength e is small enough, the motion of 
each element seems to be independent from the others. Even in such cases, the motion in 
macroscopic variables counterintuitively does not vanish in the thermodynamic limit (N — ► 
oo). This has been studied as "collective motion" in GCM |6Hl4|1, which implies some sort 
of coherence between elements. 

Fig.l gives an example of the collective motion in GCM flU) without the external 
noise (a = 0). We adopt the mean field, 

1 N 

J»»=Tr £/(*»(<)), (2) 
iV i=l 

as a macroscopic observable. While the microscopic motion shows high dimensional chaos 
in the sense that the Lyapunov dimension is proportional to the number of elements N, 
the macroscopic motion shows a quasiperiodic-like structure as is shown in Fig.l. In al- 
most all the parameter values, the mean field motion shows some coherent structure rang- 
ing from quasiperiodic-like to higher dimensional one distinguishable from random mo- 
tions ||12|| . However, even if the macroscopic motion looks like quasiperiodic, scattered 
points around the torus-like structure depicted in Fig.l does not vanish even in the thermo- 



dynamic limit |13|Jl^| , suggesting high dimensionality of the collective motion. So far, the 
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mean field motion is considered to be an infinite dimensional motion even when the torus 



like structure is observed 10-13 



The addition of noise may, however, destroy such coherence among elements. One of the 
authors reported that the microscopic external noise leads the mean square deviation (MSD) 
of the mean field distribution decreases with N ||. (The MSD decreases in proportion to 
1/N 13 with < 1, when a is larger than a certain constant.) On the other hand, it is also 
reported that, the external noise sharpens the peak in the power spectrum of the collective 
motion ||. In this Letter, we clarify the effect of noise on the collective motion in GCM. 

Consider a one-body distribution function p n (x) of the elements to study the behavior 
of the collective motion in the thermodynamic limit N — > oo. Since the mean field value 

h n = J f(x)p n (x)dx, (3) 

is applied commonly for each element, and since the additive noise can be represented as a 
deterministic diffusion process of the distribution function in the thermodynamic limit, the 
evolution of p n (x) obeys Perron- Frobenius equation written as 

f 1 (Fn(«)-z) 2 

Pn+i{x) = / dy e S7> p n (y), (4) 
J \/27T<j 



with F n (x) = (1 - e)f(x) + eh n . 

Fig. 2 gives an example of return map of the mean field value obtained numerically in 
GCM with the external noise. The parameters a and e are the same as in Fig.l. Nu- 
merical calculation was carried out through integration of Eq.(fD using a sufficiently large 
dimensional vector to approximate p n (x). As is shown in Fig.2(a), the motion has a clearer 
structure than the motion without the noise. By increasing a, motions on a torus, locking 
states and lower dimensional chaos are observed (Fig.2(b)). With the further increase of a, 
the collective motion collapses to a fixed point. Hence, with the increase of the noise a sort 
of bifurcation to lower dimensional motions is observed. 

To clarify this point, it seems natural to measure the number of effective degrees of 
freedom in the collective motion. We calculate the Lyapunov dimension of the dynamics of 
p n (x). The Lyapunov exponents are given by the growth rates of tangent vectors around 
the orbit of Eq.(f|). For numerical calculation, p n (x) is approximated by a sufficiently large 
dimensional vector, and its linear stability around the orbit is studied. 

In Fig. 3, the Lyapunov dimension denoted by D c is plotted as a function of the noise 
variance a 2 . For sufficiently large a, only the stationary state is observed and Dq is zero 
accordingly. With the decrease of a we have found the low dimensional collective mo- 
tion (Dc ~ 0(1)) ranging from the motion on a torus to low dimensional chaos. With the 
further decrease of a, the dimension grows as 

D c oc -log a 2 . (5) 

This implies that the number of effective degrees of freedom goes to infinity in the zero noise 
limit, as is expected from the analysis of the collective motion in GCM without the external 
noise. 

In the large a regime a variety of bifurcations appears, which may strongly depend on 
the parameters. However, the above result suggests that the scaling relation (|5|) will be a 
characteristic common to the high-dimensional collective motions in the small a regime. 
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Although the evolution rule is originally given for the microscopic variables, our main 
interest is on the behavior of macroscopic variables which would be the only possible ob- 
servable in typical cases. Thus, it is highly desirable to obtain a closed description of the 
behavior of the macroscopic variables, which could be written as 

h n = h(h n -i,h n -2,---), ( 6 ) 

for an idealized example. In most cases, however, it is quite difficult and may well be 
impossible to obtain such a description. Thus, we examine the linear response of the system 
against infinitesimal perturbation on the macroscopic variables, and obtain the variational 
equation describing the evolution of the small deviation of the macroscopic variables in a 
neighborhood of a trajectory. 

In the present case, since the elements interact only through the mean field value, it is 
quite natural to expect that the behavior of the mean field value can be consistently described 
by itself. We expect that the effective number of the dimension of the collective motion gives 
substantial agreement with the Lyapunov dimension of the macroscopic dynamics estimated 
in the above mentioned way. We concentrate on the small a regime and give qualitative 
explanation for the scaling relation (|5|). 

If we consider small deviations r] n of h n , then r] n is regarded as a function of 
{Vn-i, Vn-2, ■ ■ ■}• The evolution of i] n from the unperturbed orbit h n is given by, 

oo 

r] n = J2 L rVn~T + 0(r] 2 ), (7) 

T=l 

where L T is a coefficient to give the linear response of the mean field value at n step to 
the displacement at n — r step. The number of the Lyapunov dimension of the mean field 
dynamics is estimated from the eigenvalues of this linear regression |16| . 



First we estimate L T from the dynamics of the distribution function given by Eq.(|j). In 
the small noise limit (a — > 0), from Eqs.(^) and (|j), L T is given by 

rf^ft^W, (8) 



dx 



where F«(x) e^o^o-o F„_ 

For r ^> 1, dF£'(x)/dx in Eq.(§) changes its sign quite frequently in x. Now let us 
consider the partition of x at the points such that F^{x) = 0. Denoting the typical value 
of dF^ T \x)/dx by d(r), the interval of partitions is estimated at l/d(r), which decreases 
rapidly with r. Since the integration in a partition becomes zero if p(x) stays constant in that 
partition, the partitions where p n (x) changes drastically in x contribute to the estimation 
of L T much more than the partitions where p n (x) does not changes so much. 

In the case of small noise limit (a — > 0), the most drastic change of p n (x) comes from the 
inverse square-root singularities, which is the characteristic structure of distribution function 
for the logistic map. Hence, the integration in the partitions containing the characteristic 



structure in p n (x) is estimated as 0(Jd(r)) |0|. d(r) is roughly estimated at e XmT for 



r ^> 1, where A m is the Lyapunov exponent of the local mapping [18]. Consequently, the 



response L T to the perturbation grows exponentially with the rate \\ r 
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Even in the presence of finite amplitude of the noise, the above order estimation for L T 
is still valid for r smaller than r c = — logcr/A m , where the typical width of the partitions 
becomes comparable with the typical amplitude of the noise, i.e. l/d{r) ~ e _AmT = a. 

For larger r > r c , however, the effect of noise in smoothening the distribution p(x) 
appears so that L T will start to decay with r. 

Partially integrating (ffi), we obtain 



dxF^\x) dpn - T + l{x) 



1-eJ n v ' dx 

dxF^\x) dpn -]^ X \ (9) 



1-eJ n v ' dx 

with F^\x) = F^\x) — F^ , where F^ is the average value of F^\x) over the support 
of Pn(x). With the increase of r, F^\x) becomes a rapidly oscillating function about zero- 
mean in x, and the integration of F^(x) over any finite range within the support of p n (x) 
will approach zero [|17|]. Since | dpn -[+^ x ) | j s uniformly bounded due to the existence of the 



noise, L T converges to zero. Hence, as far as the linear stability is concerned, the mean field 
value is not sensitive to the mean field values before sufficiently long time step (before 0{r c ) 
steps). Thus we can consider a dynamics of h n as a function of the mean field values of the 
past 0(t c ) steps as Eq.(^) at least in the neighborhood of the orbit. 

In summary, the amplitude of L T grows with r as L T ~ C 2 AmT for r < 0(r c ) whereas 
it starts to decay for r > 0(r c ). Thus for sufficiently large r c , i.e., for sufficiently small a, 
the number of positive eigenvalues around zero for Eq.(|7|) is estimated at 0(r c ). Since we 
have to consider the contribution only from the latest 0(t c ) steps, the dimension of this 
dynamical system can be at most 0(r c ). Hence, the dimension of the mean field dynamics is 
within the order of r c . Accordingly the number of effective degrees of freedom of the mean 
field dynamics grows as — log cr with the decrease of a, and can grow arbitrary large as a 
approaches zero. 

In the present Letter, we have shown that the noise in a microscopic level reduces the 
complexity of the collective motion, which is characterized by the number of degrees of 
freedom. It is shown numerically that the dimension of the dynamics of Perron-Frobenius 
equation Eq.(f|) satisfies On the other hand, analysis on the mean field dynamics also 
supports ([)]). Hence, the number of effective degrees of freedom of the collective motion in 
the present GCM is concluded to satisfy the scaling relation (|5|). 

Such a relation is expected to hold when the collective motion keeps high dimensional 
motion with microscopic chaos. High-dimensional collective motions are supported by the 
exponential growth of the linear response coefficient. When the distribution has peak con- 
centrations within a small width 6x, such an exponential growth is expected until r ~ log5x 
and thus a high dimensional collective motion appears f20fl . Even when the local dynamics 
is given as a higher- dimensional non-hyperbolic map or a chaotic flow, the present argu- 
ment on the collective motion is expected to hold. Hence, the logarithmic dependence of 
the dimension upon the noise intensity given by (|5|) will be observed in a broad range of 
systems. 

The induced regularity by the addition of noise was also reported as Noise-Induced 
Order (NIO) in a one-dimensional map 0. The mechanism to induce regularity in our system 
has a similarity with NIO case, in the point that the external noise destroys a dynamical 
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structure which causes the irregular behavior. In NIO, the noise reduces the concentration of 
measure in the instability region where the intermittent behavior is generated. On the other 
hand, in the present case, the noise smoothes the singularity of the distribution function, 
that is the source of high- dimensional instabilities of the collective motion. 

In the case of NIO, however, the motion is not completely regular, since the noise is 
imposed upon the observed variable itself. On the other hand, in the present case, the 
noise is imposed on the microscopic variables, whereas observed is the macroscopic collec- 
tive variable. Thus, the noise-induced motion is deterministic and low dimensional in the 
thermodynamic limit. 

We should also mention that our present result may be applicable to experimental sys- 
tems, such as fluid turbulence and neural systems consisting of a huge number of neurons 
with nonlinear dynamics. By controlling the thermal noise or some external fluctuations, we 
expect to observe the gradual reduction of the dimension and noise-induced low- dimensional 
motions, and hence the existence of the collective chaos will be clarified. 
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FIGURES 

FIG. 1. A return map of the collective motion in GCM (|l]) without noise. 
a = 1.86, e = 0.1, a = 0.0, N = 10 7 . 



FIG. 2. A return maps of the collective motion in GCM (|l]) with noise, a = 1.86, e = 0.1. 
The noise variances are (a) a 2 = 1.5 x 10~ 6 , and (b) a 2 = 3.0 x 10~ 6 . Numerical calculation was 
carried out with integration of Eq. (|j) using a sufficiently large dimensional vector to approximate 
the distribution function. 



FIG. 3. The dimension of the collective motion is plotted as a function of noise intensity, ob- 
tained by the Lyapunov dimension. The Lyapunov exponents are calculated as growth rates of the 
tangent vector around the orbit, obtained by Eq.(|j). For numerical calculation, p n {x) is approxi- 
mated by a sufficiently large dimensional vector, and the tangent vectors are orthonormalized at 
each time step. 
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